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Abstract 

The thrust augmentation of a set of ejectors driven by a shrouded Hartmann- Sprenger tube has been measured at 
four different frequencies. Each frequency corresponded to a different length to diameter ratio of the pulse of air 
leaving the driver shroud. Two of the frequencies had length to diameter ratios below the formation number, and 
two above. The formation number is the value of length to diameter ratio below which the pulse converts to a vortex 
ring only, and above which the pulse becomes a vortex ring plus a trailing jet. 

A three level, three parameter Box-Behnken statistical design of experiment scheme was performed at each 
frequency, measuring the thrust augmentation generated by the appropriate ejectors from the set. The three 
parameters were ejector length, radius, and inlet radius. 

The results showed that there is an optimum ejector radius and length at each frequency. Using a polynomial fit 
to the data, the results were interpolated to different ejector radii and pulse length to diameter ratios. This showed 
that a peak in thrust augmentation occurs when the pulse length to diameter ratio equals the formation number, and 
that the optimum ejector radius is 0.87 times the sum of the vortex ring radius and the core radius. 


Nomenclature 

a radius of vortex ring core 

A p ui se exit area of the driving jet 

Dpuise ^A pu i se /n (44 mm) 

/ frequency of driving jet 

K circulation in vortex ring 

L e j length of the ejector 

L pu i se length of the slug of air in each pulse 

riije, mass flow rate of the jet 

nipuise mass of air in each pulse 

m tota i total mass flow rate leaving the ejector 

N formation number 

p fraction of total jet momentum in the vortex ring 

R radius of the vortex ring 

R e j radius of the ejector 

R ex it exit radius of the ejector 

R s inner radius of the shroud (25.4 mm) 

r„ ejector inlet radius 

S t Strouhal Number (f D pu i se /U) 

t time 

t ’ time at which the pulse ratio L puls JD pu i se =N 

Tjet thrust of jet alone 

u ( t ) jet velocity as a function of time 

U average velocity of the driving jet = Tj et lrhj et 

Vexit average in time and space of ejector exit velocity 

a thrust augmentation coefficient 

a ss thrust augmentation for L pu i se /D pu i se = oo 
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o Uotai weighted sum of a ss and a,,,- 

a vr thrust augmentation due to vortex ring only 

a v i- y a vr for Lp U i se l D p U i se N 

P mass entrainment coefficient = {m tota i - rhj et )/ rhj et 

e a/R 

^ M jet /P U ^ pulse 

p atmospheric air density 

t period of one pulse = 1 If 


I. Introduction 

The current interest in Pulsed Detonation Engines has revived interest in pulsed ejectors (refs. 1 to 3), since the 
latter can increase the thrust, and potentially reduce the noise level, of such engines with little additional hardware. 
There have been measurements on pulsed ejector performance in the past, using drivers that were not detonation 
devices. Lockwood (ref. 4), using a pulsejet driver, measured thrust augmentations as high as 1.9, with a very short 
ejector. Didelle (ref. 5), using a jet interrupted by a butterfly valve, also observed a maximum thrust augmentation 
ratio of 1.9, but with a much longer ejector. Bertin (ref. 6), also claimed a thrust augmentation ratio of around 2, but 
gave few details. More recently, Paxson, Wilson, and Dougherty (ref. 7) have repeated Lockwood’s experiment, 
with similar results, although the optimum ejector was rather longer than was Lockwood’s. Wilson and Paxson 
(ref. 8), using a Hartmann- Sprenger, or resonance, tube to generate a pulsed jet, measured a maximum thrust 
augmentation of 1.32, and found, as did Paxson, Wilson and Dougherty, that there was an optimum radius for the 
ejector. Despite these efforts, it is still not clear what are the important parameters for optimizing a pulsed ejector. 
Except for that of Didelle, in none of the above experiments was the length of the pulse driving the ejector changed. 
Didelle, who varied the jet frequency, and hence pulse length, found his maximum thrust augmentation at a Strouhal 
number of 0.07. 

A major feature of any pulsed jet flow is the starting vortex ring (Elder and deHaas (ref. 9), Das, Arekeri, 
Krothapalli, and Lourenco (ref. 10)). Gharib, Rambod, and Shariff (ref. 11) have shown that there is a critical value 
of the pulse length to diameter ratio L pu i s JD pu i se , which they called the formation number N: for values of 
LpuisJDpuise < N, the pulse transforms into a vortex ring only, for L pu i se /D pu i se > N, the flow becomes a vortex ring 
followed by a trailing jet. Since the structure of the flow changes in this way with L pu i se ID pll i se , one would expect that 
changing L pu iJD pu i se might also affect ejector performance. 

In order to explore this possibility, experiments to measure pulsed thrust augmentation were made using a 
shrouded Hartmann-Sprenger tube as the pulsed source. This source, previously used by Wilson and Paxson (ref. 8), 
has the advantage that by varying the tube and shroud lengths, different frequencies can be generated, corresponding 
to different values of the pulse length to diameter ratio. A statistical experiment was designed to explore the effects 
of ejector length, radius, and nose radius. This was repeated at each of four frequencies, and by fitting the data with 
a two parameter model, the effect of L pu i sl JD p „t se and ejector radius on the thrust augmentation was determined. The 
maximum thrust augmentation was found to be at an L pu i se /D pu i se equal to the formation number, and the optimum 
ejector radius was found to be slightly smaller than the size of the vortex ring, defined as ( R + a). 

II. Apparatus 

A photograph of the apparatus is shown in figure 1. Consisting of a shrouded Hartmann-Sprenger tube, a thrust 
plate, and a set of ejectors, it is the same as described by Wilson and Paxson (ref. 8), where a more detailed 
description will be found. The Hartmann-Sprenger tube (fig. 2) is a device which converts a steady, supersonic flow 
into a pulsed flow, by directing it into a resonance tube. It was shrouded to channel the pulsed flow into the axial 
direction. As initially built, the Hartmann-Sprenger tube operated at a frequency of 550 Hz. Different frequencies 
were obtained by, first, putting a plug inside the resonance tube to shorten it, which produced a frequency of 
1100 Hz, and then making two longer tubes and shrouds to replace the original tube, to give frequencies of 275 and 
125 Hz. Since the nozzle upstream of the resonance tube is choked, the mass flow remained the same for all 
frequencies. 

The set of ejectors is also given in figure 2, to illustrate the geometry. The set includes four different ejector 
radii, R cj , at each of which three lengths, L (:j , are possible, as well as three inlet radii, r n . Since there was a diffuser on 
each ejector, the exit radius is greater than the ejector radius itself. 
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In running the experiment, the thrust of the jet without an ejector was measured, followed by the thrust with an 
ejector. The ratio of the two thrusts is the thrust augmentation a. The exact experimental procedure is described in 
Wilson and Paxson (ref. 8). 



Figure 1. Photograph of the apparatus, showing the thrust plate at the top, the shrouded 
Hartmann-Sprenger tube at the bottom, and an ejector above it. 
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Pulsed flow out 



♦ 

Steady flow in 


Figure 2. Schematic diagram of the Hartmann-Sprenger tube 
and the set of ejectors used in the experiment. 


ITT. Design of the experiment 

The original experimental design of Wilson and Paxson (ref. 8) for a frequency of 550 Hz was chosen to be a 
statistical design of experiments, three level, three parameter Box-Behnken scheme (ref. 12), with L e j/R s , R e jlR s ,, and 
r,JR s as the three parameters. This does not mean to imply that division by R s is an appropriate scaling; it is simply a 
convenient way of making the numbers non-dimensional. Initially, the values of R ej /R s were 1.1, 1.5, and 3.0. As the 
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thrust augmentations for R ej /R s =3.0 were small, it was found useful to add an additional group of ejectors with 
R e jlR s — 2.0, using R e j/R s = 1.5 as the “central” group. In making runs at a frequency of 1100 Hz, it was found that 
the optimum diameter for thrust augmentation appeared to be shifting towards lower values of R e j/R s , and so runs 
were made for R e j/R s = 1.1, 1.5, and 2.0 only, and not 3.0. As the frequency was lowered, the optimum shifted 
towards higher R e j/R s , and so the values of R e /R s used were 1.5, 2, and 3, and not 1.1. The full set of runs made is 
given in table 1 . 


TABLE 1.— THE MATRIX OF TEST RUNS COMPRISING THE 3 PARAMETER, 

3 LEVEL STATISTICAL EXPERIMENT AT EACH OF THE FOUR FREQUENCIES USED 


Run 

Rej/Rs 

Lej/Rs 

rJRs 



a 

No. 




/= 1100 

/= 550 

/= 275 /= 125 

1 

1.1 

3.125 

0.5 

1.17 

1.155 


2 

1.1 

7.125 

0.25 

1.149 

1.16 


3 

1.1 

7.625 

0.75 

1.134 

1.155 


4 

1.1 

12.375 

0.5 

1.122 

1.092 


5 

1.5 

2.875 

0.25 

1.099 

1.264 


6 

1.5 

3.375 

0.75 

1.053 

1.275 


7 

1.5 

7.375 

0.5 

1.172 

1.330 


8 

1.5 

7.375 

0.5 

1.167 

1.325 


9 

1.5 

7.375 

0.5 

1.169 

1.308 


10 

1.5 

12.125 

0.25 

1.166 

1.255 


11 

1.5 

12.625 

0.75 

1.165 

1.267 


12 

2 

3.125 

0.5 

1.022 

1.118 


13 

2 

7.125 

0.25 

1.092 

1.206 


14 

2 

7.625 

0.75 

1.10 

1.226 


15 

2 

12.375 

0.5 

1.167 

1.266 


16 

3 

3.125 

0.5 


1.034 


17 

3 

7.125 

0.25 


1.075 


18 

3 

7.625 

0.75 


1.085 


19 

3 

12.375 

0.5 


1.139 


20 

1.5 

3.125 

0.5 



1.287 1.14 

21 

1.5 

7.125 

0.25 



1.320 1.221 

22 

1.5 

7.625 

0.75 



1.355 1.237 

23 

1.5 

12.375 

0.5 



1.326 1.265 

24 

2 

2.875 

0.25 



1.277 1.186 

25 

2 

3.375 

0.75 



1.343 1.218 

26 

2 

7.375 

0.5 



1.379 1.317 

27 

2 

7.375 

0.5 



1.373 1.295 

28 

2 

7.375 

0.5 



1.376 1.313 

29 

2 

12.125 

0.25 



1.305 1.271 

30 

2 

12.625 

0.75 



1.383 1.310 

31 

3 

3.125 

0.5 



1.053 1.075 

32 

3 

7.125 

0.25 



1.106 1.137 

33 

3 

7.625 

0.75 



1.116 1.182 

34 

3 

12.375 

0.5 



1.197 1.250 


IV. Experimental results 

The measurements of thrust augmentation found in the various sets of Box-Behnken runs are given in table 1 and 
in figures 3 (a) to (d), in which thrust augmentation is plotted against normalized ejector length, LJR S , for each 
normalized ejector diameter, R e j/R s , at different values of the normalized nose radius, rJR s . The 90 percent 
confidence level in a is ± 0.03. The data from a Box-Behnken three parameter set can be fitted with a response 
surface of the form 
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Figure 3. Results from the statistical experiments at four different frequencies, (a) 1100 Hz. (b) 550 Hz. 

(c) 275 Hz. (d) 125 Hz. Different symbols are used for each nose radius, ■, r n /R s = 0.25; ♦, r n /R s = 0.50; 

* , r n /R s = 0.75. Different colors are used for each ejector radius, red, R e j/R s = 1.1; dark blue, R e j/R s = 1.5; 
orange, R e j/Rs = 2.0; green, R e j/R s = 3.0. The points are experimental data, the lines are response surface 
fits: dotted, r n /R s = 0.25; solid, r n /R s = 0.50; dashed, r„/R s = 0.75. 


+ b x ■ ( L ej /R s ) 

+ b 2 ■ ( R e j /R s ) 

+ b 3 ■{>'„/ Rs ) 


+ i • ( L e j / R s ) 2 

+ b 22 ■ {R e j /R s ) 2 

+ ^33 •( r n/R s )~ 

(1) 

+ - (b e j /R s )• {r c j /R s ) 

+ ^13 ' {b e j /R s )■ ( r n /Rs) 

+ b 22 - ( R e j /R s ) ■ (r n /R s ) 



where the values of the constants b,y are determined from the data. These were found by inserting the data into a 
computer program (Seshadri and Demming (ref. 13)), which calculates the constants, and values of the confidence 
level for each constant. Constants with low confidence level were eliminated, until only terms with levels greater 
than 90 percent were retained. The resulting constants are listed in table 2. Once the response surface is known, the 
optimum thrust augmentation can be predicted, together with the values of L e j/R s , R e j/R s , and r n /R s at which the 
optimum occurs. These values are also listed in table 2. 
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TABLE 2.— CONSTANTS FOR EQUATION (1)— THE RESPONSE SURFACE OF THE BOX-BEHNKEN 
EXPERIMENTAL DESIGN, TOGETHER WITH PREDICTIONS FOR THE 
MAXIMUM THRUST AUGMENTATION AT EACH FREQUENCY. 



/= 1100 Hz 

/= 550 Hz 

/= 275 Hz 

/= 125 Hz 

b 0 

1.16812 

-0.18464 

0.54744 

0.32863 

bi 

-0.01464 

-0.00713 

- 

0.03207 

b 2 

— 

1.81433 

0.68699 

0.64294 

b 3 

0.28836 

0.23637 

0.52521 

0.44842 

bj i 

-0.00126 

-0.00197 

-0.00103 

-0.00193 

b 22 

-0.07511 

-0.62527 

-0.19356 

-0.15947 

b 33 

-0.41306 

-0.21744 

-0.30355 

-0.39603 

bi 2 

0.02336 

0.02591 

0.01052 

0.00505 

bi 3 

0.01187 

- 

- 


b 23 

- 

- 

-0.06310 


^max 

1.189 

1.331 

1.385 

1.338 

Rej/Rs 

1.98 

1.62 

1.91 

2.21 

L e j/R s for C^max 

12.6 

8.9 

9.9 

11.2 

r n /R s 

0.5 

0.5 

0.75 

0.5 

Le/Rei 

6.36 

5.49 

5.18 

5.07 


Although equation (1) includes r„, changes in a due to changes in r„ for the range of values used are very small, 
and the changes seen experimentally are statistically insignificant, except at a frequency of 275 Hz where for 
R e /R s = 2, thrust for r,JR s = 0.25 seems lower than values for r n IR s = 0.5, and r,JR s = 0.75 indicating an optimum 
value between the latter two values. That, apart from this one case, changes in the nose radius seem unimportant is 
surprising, and may be a consequence of the limited values of r n used in this experiment. Other experiments (refs. 7 
and 15) have shown a significant effect of r n , including the possibility of an optimum value. Sections through the 
response surface at each value of ejector radius, R e j, for rJR s = 0.5, are also plotted in figure 3, showing good 
agreement with the experimental results. The optimum values predicted are little different from the maximum 
observed experimental values. 

V, Interpretation of the experimental results 
A. Statement of the problem 

It was originally hoped that by finding a response surface at each frequency, some insight would be gained into 
how the parameters affect the thrust augmentation as a function of the frequency. Except for the conclusion that 
changes in nose radius were not statistically significant, this was not apparent from the constants listed in table 2. 
Since the objective here is to find the overall maximum thrust augmentation, and at what values of the parameters 
this occurs, only the maximum augmentation achieved at any diameter, length and frequency will be considered 
further. These are listed in table 3. A more detailed look at these data is required. 

TABLE 3 .—L pu iJD pulse AND s FOR EACH FREQUENCY USED, TOGETHER WITH THE MAXIMUM THRUST 
AUGMENTATION OBSERVED FOR EACH R ej , TOGETHER WITH VALUES OF THRUST AUGMENTATION 
FOR L pu i se /D pu i se = N, I.E., a vr |,v, AND FOR L pll iJD puhe = co, I.E., a ss , AS CALCULATED FROM 
THE VALUES OF a AT /= 275, AND/=125. 


Frequency (Hz) 

1100 

550 

275 

125 


0 

LpulsJDpuise 

2.5 

5 

10 

22 

6.86 =N 

00 

8 

0.177 

0.386 

0.54 

0.54 

0.54 

0.54 

R/R s 

1.2438 

1.313 

1.3084 

1.3084 

1.3084 

1.3084 

max for ReJR-s — 1 • 1 

1.17 

1.16 

- 

- 

- 

- 

W'max for ReJR-s — 1*5 

1.169 

1.321 

1.355 

1.265 

1.395 

1.191 

^ max for ReJR-s 2.0 

1.167 

1.266 

1.376 

1.308 

1.406 

1.252 

<W for RJR S = 3.0 

- 

1.139 

1.197 

1.25 

1.174 

1.294 


B. Relationship of frequency to L plUse ID pu i se 

The experiments were performed at four different frequencies. Since the mass flow was constant at all 
frequencies, the mass per pulse m pu i se is given by 
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M pulse W jet / f 

But m p is also given by the product of pulse volume and density 


(2) 


M pulse P tt D p U i se ^ pulse /4 (2) 

from which it is seen that L pu i se /D pu i se is inversely proportional to frequency. In table 3 the values of L pu ) se /D pu i se for 
each frequency used are given. Using the slug model plus some experimental measurements Wilson (ref. 14) has 
indicated that the value of the formation number is 6.86 for this shrouded Hartmann-Sprenger tube. Thus the two 
higher frequencies have values of L pu i se /D pu i se below the formation number, whilst the two lower frequencies have 
values of L pu i s JD pu j se above it. 

C. Thrust augmentation for L pu i se ID pu i se greater than N. 

For the measurements at values of L pu i se /D pu i se above the formation number, the flow consists of a leading vortex 
ring and a trailing jet. The vortex ring has exactly the same size and properties as that generated at L pu i se /D pu i se = N. 
The thrust on the ejector will be partially due to the vortex ring and partially to the trailing jet. There will be a time, 
t ' , following initiation of the jet at which the value of L pu i se = \u(t)<lt is equal to N.D pu i se . The fraction of the total jet 
impulse that is in this portion of the jet is 


P = 



(4) 


It is this fraction of the impulse that goes into forming the vortex ring, giving rise to a thrust augmentation a vr . The 
remaining impulse will be in the trailing jet, which will be assumed to create a thrust augmentation a ss as would a 
steady jet. Thus, for L pu i se /D pu i se greater than the formation number, the total thrust augmentation will be 

®-total = p CL w + (l — ss (5) 

in which a w is calculated at the formation number, i.e., it is a, r |.\/. A measurement of the jet exit velocity u(t) for this 
tube has been reported previously (Wilson (ref. 14)). Given this distribution, the integrals which determine p can be 
performed numerically. The result is given in table 4. Of course, if L pu i se /D pu i se < N, then a tota i = a,,,, with no 
contribution from a ss . 


TABLE 4.— VALUES OF p AS A FUNCTION OF L pu JD pulse 

FOR THE JET VELOCITY OF THE HARTMANN-SPRENGER TUBE 

L pu iJD pulse 6.86 7.4 8 10 12 14 16 18 20 22 24 oo 

p= 1 0.9736 0.9424 0.8064 0.6724 0.5742 0.5002 0.4443 0.3990 0.3649 0.3363 0 


At each of the two low frequencies, the value of p is known from table 4, i.e., p = 0.8064 at L pu iJD pu i se = 10, and 
0.3649 at L pu i se /D pu i se = 22. At each frequency, there are three values of a tota u one for each different ejector radius. 
Thus for each radius, there are two values of a tota h which can be substituted into an equation (5) for each frequency, 
so that the resulting pair of equations can be solved for a„.|AS and a. s . s . Values of a w .| 7 v, and a ss are given in table 3. 
The values of a ss for each ejector radius are fitted by the equation: 

a ss = 1.094 + 0.158 In (2 R exit / D pulse ) (6) 

As can be seen, a ss increases as the ratio of ejector exit diameter to tube diameter increases. This trend, and the 
values of the thrust augmentation are in agreement with values of steady state thrust augmentation given by Porter 
and Squyers (ref. 16). 

D. Thrust augmentation for L pu i s JD pu i se less than N 

The objective of this work is to gain some understanding of the conditions which lead to maximum thrust 
augmentation, in particular ejector radius, and vortex ring size. Length is of course also an important parameter, but 
the value of L e j ecto JR e j ector at a max , also given in table 2, and plotted in figure 4, does not appear to be very dependent 
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on LpuiJD pu i se . At L pu iJD pldse = N, the value of L ejector IR eJector for a max is 5.2 - 5.5, and since from figure 3, at the 
maximum value of thrust augmentation for any of the curves, the augmentation does not depend strongly on ejector 
length, a value of L e j ec toJRejector = 5.3 would probably be reasonably close to optimum at any frequency. To 
investigate the effects of ejector radius and vortex ring size, it was decided to use the values of maximum thrust 
augmentation observed at the two higher frequencies, together with the values calculated for o, w \n from the low 
frequency data, as the dependant parameter in creating a response surface in parameter space. A major question was 
what are the appropriate independent parameters to use in forming a response surface. Clearly, the thrust 
augmentation depends on the vortex ring properties, and the ejector size. Although L pu i se /D pu i se is important for 
vortex ring formation, once the ring is formed, the ratio e is what characterizes it, and was chosen as the parameter 
for the vortex ring. For the ejector, the choice was not so clear. Considering the thrust augmentation, if the total exit 
flow from the ejector can be described by an average velocity in time and radius, V CXII , and the jet has an average 
velocity U such that Tj et = m jet U, then: 

a = total thrust/thrust of jet alone 

— total Vexit jet U 

= pn{R exi tV ex it) 2 /pnt,(D pu i se U/lf 

= (2R exit /D pulse f {V exit /U) 2 k 

from which: 

Kxit /U = Ja£,{D p U is e /2 R exit ) 


( 7 ) 


( 8 ) 
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Dpulse^2R e xit 

Figure 5. V cx jj/U versus I ^pulsi/ ~ R cxi t • ® Ppulsc^RpuIsc = 2.5, ® ^ pulsi^pulsi' — 5, 

□ Lpuise/Dpuise = 6.86, calculated from measurements of a, x Lp U | se /Dp U | se = 5, 
from direct measurement of V ex j t /U. 


Values of V ex j,/U derived from equation (8) are plotted against D pu i s J 2R ex i t for the points in table 3 with 
LpuiselDp U i se less than or equal to N in figure 5, as well as two measured values of V exil /U reported by Wilson and 
Paxson (ref. 8). Within the experimental error, these points are all fitted by a single line, shown in the figure. The 
equation of the line is: 

Vexit/U = -0.123 + l.442D pulse /2R exit - 0M%{D pulse /2R exll ) 2 (9) 

This suggests that D pu i s J2R exit would be an appropriate parameter for fitting V ex i t IU to the ejector size, from which 
the thrust augmentation could be obtained. Unfortunately, this did not prove to be the case, since small errors in 
fitting V ex iJU are magnified in calculating thrust augmentation, as it is the square of V ex iJU that appears in the 
equation for thrust augmentation. Nevertheless, the fit to V ex i t IU is useful in calculating an estimate of the mass 
entrainment, since, as can be shown easily: 

a = (l + P) Vexit/U (10) 

Instead of using the fit to V ex JU, all the values of thrust augmentation given in table 3 were entered into a non- 
linear regression analysis with s and one of either D pu i s J2R exih RIR ex u, or (R+a)/R lx i, as the independent parameters. 
It was found that the best fit was with (R+a)/R ex j t as the second parameter, with the resulting fit being: 

a = 1.816-3.743s + 0.0389s 2 - 1.447 (R + a)/R exit + 0.679 ((R + a)/R exit f ^ ^ 

+ 8.43 Is (R + a)/R exlt - 4.014s {{R + a)/ R exi , f 

Using this fit, together with thrust augmentation calculated from equation (5) for L pu i se /D pu i se greater than N, 
thrust augmentation is plotted against ejector radius in figure 6, and against L pu iJD pu i se in figure 7. In using this fit. 
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the values of R and e (and hence a) as a function of L pu i se ID pu i se were derived from slug theory as described by 
Wilson (ref. 14), R /R s and s are given in table 3. 

VI. Discussion 

The plot of thrust augmentation versus ejector radius in figure 6 shows that the optimum ejector radius increases 
as L P uisJ D pu i se increases, as had already been observed in figure 2. The optimum ejector radius for L pl ,i se /D pu i se = N is 
at R e j/R s = 1.75, which corresponds to R e J{R + a) = 0.87, i.e., the optimum ejector is slightly smaller than the vortex 
ring size. The plot of thrust augmentation against L pu j se !D pll i se in figure 7 shows that there is a peak thrust 
augmentation at L pu i s JD pu i se = N. However, as L pu i se /D pu i se gets very much larger than N, the thrust augmentation 
approaches a„, which increases monotonically with R ex u, so there will no longer be an optimum ejector radius at 
very large L pu iJD pu i se . The optimum thrust augmentation at L pu i se /D pu i se = N is around 1.4, considerably lower than 
the values around 2.0 reported by Lockwood (ref. 4), Didelle (ref. 5), and Paxson, Wilson, and Dougherty (ref. 7). It 
has been shown (Wilson (ref. 14)) that the circulation produced by the annular exit of the shrouded Hartmann- 
Sprenger tube ( K = 13.7m 2 /s) is lower than would be expected for the same value of U from a circular exit ( K = 20 
m 2 /s). Since the vortex ring velocity is proportional to the circulation, a circular exit would have a higher value of 
vortex ring velocity, and hence, presumably, also a higher value of V exiL If it is also assumed that the thrust 
augmentation is proportional to V ex n IU, then a circular tube producing the same exit velocity as the Hartmann- 
Sprenger tube coidd be expected to generate a thrust augmentation of 2.04, in line with the values from the other 
experiments. However, equation (7) shows thrust augmentation proportional to (V exit /U) 2 , which woidd suggest that 
a maximum a,,,- of 3 should be possible. 

These experiments indicate that a length of about two and a half times the ejector diameter is optimum. This is 
longer than claimed by Lockwood (ref. 4), but shorter than optimum length to diameter ratios found by Paxson, 
Wilson and Dougherty (ref. 7), or Didelle (ref. 5), The results given in table 2 for the optimum ejector length to 
diameter ratio seem to show a gradual decrease in the ratio as L pu \JD pu i se increases. Thus the discrepancy in 
optimum ejector length to diameter between this work and those of Lockwood, and Paxson, Wilson and Dougherty 
might be explained as due to being in very different regimes of L pu i se ID pu i se . Didelle’s maximum thrust augmentation 
was at a Strouhal number of 0.07. Because his valve gave two pulses per second, the value of L pu iJD pu i se used by 
Didelle is 1/2 S t , i.e., 7 at maximum thrust augmentation. Most of his runs were for L P uisJD pu is e < 7, with a few at 
LpuiselD pu i se = 10. It is gratifying that Didelle’s optimum value of L pu i se /D pll i se is in close agreement with that found 
here, but that leaves the discrepancy of the optimum ejector length to diameter ratio unexplained. 
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Pulse length/diameter Lp U | Se /Dp U | Se 


Figure 7. Thrust augmentation versus Lp„i se /Dp U ] se for each value of R e j/R s . The 
thin horizontal lines at the right represent the values of a ss . The data and fitted 

lines correspond to different values of R e j/R s , namely: R e j/R s = 1.1, — , *, 

Rej/Rs = 1*5? ' ' Rej/Rs = ^ Rej/Rs = 3, , D. 


In these experiments, Didelle’s experiments, and the experiments with pulsejet drivers, the length of the pulse 
increases with the inverse of the frequency, and so low frequency approaches steady state flow. Both Didelle 
(ref. 5), and Wilson and Paxson (ref. 8) found that V exit had a large average component, with an oscillation on top of 
it, such that the ejector flow never reached zero velocity. Care must be used in relating the results to pulse 
detonation devices, where the pulse length is determined by the length of the tube, and is independent of frequency. 
Because the duty cycle is very low in detonation devices, V ex n may drop to zero between pulses. This is less likely as 
the frequency increases, and so high frequency is more likely to approach a steady state flow in the ejector. 

VII. Conclusions 

In conclusion, for the Hartmann- Sprenger tube investigated, unsteady thrust augmentation is shown to have a 
peak when the driving pulse jet has a pulse length to diameter ratio equal to the formation number. For values of 
LpuisJDpuise close to N, the optimum ejector radius will be equal to 0.87(7? + a). As L pu i se /D pu i sc , increases, the 
optimum ejector radius also increases, until for very large values of L pu i se /D pu i se , the thrust augmentation approaches 
steady-state values, increasing monotonically with ejector radius. The optimum ejector length appears to be about 
2.5 times the ejector diameter, although this does not agree with the results of other workers (refs. 4, 5, and 7). 
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